Photovoltaic Hall effect in graphene 
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Response of electronic systems in intense lights (AC electric fields) to DC source-drain fields is 
formulated with the Floquet method. We have then applied the formalism to graphene, for which 
we show that a non-linear effect of a circularly polarized light can open a gap in the Dirac cone, 
which is predicted to lead to a photo-induced dc Hall current. This is numerically confirmed for a 
graphene ribbon attached to electrodes with the Keldysh Green's function. 
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Introduction Non-linear phenomena in electronic 
systems are fascinating since they can lead to trans- 
port properties qualitatively distinct from those in equi- 
librium. In this Letter we seek such a possibility by 
combining (i) the geometric phase argument extended to 
electron transports in intense AC fields with (ii) physics 
of graphene involving chiral states associated with two 
Dirac cones. Indeed, the geometric phase has become 
an important ingredient in the modern theory of electric 
transport [l| , which goes back to Thouless's idea of charge 
pumping where he showed that an adiabatic deformation 
of the system may lead to quantized transport 0. Ex- 
tensions along several directions have been done subse- 
quently. Thoulcss et al. (TKNN) [|[ have shown that the 
Kubo formula for the Hall conductivity can be expressed 
as a topological density. Berry then showed that such 
phases are present in general quantum systems, and the 
topological density in the TKNN formula is now called 
the Berry curvature 0]. Aharonov and Anandan further 
extended the notion of geometric phases to non-adiabatic 
situations, i.e., Aharonov- Anandan (AA) phaseQ. 

With these as a background, the concept of the present 
study, with graphene in mind, is simple. Let us put a 
crystal with a Dirac band in a circularly polarized light. 
As in Fig. [IJa), an intense AC field A ac (t) will deform 
the single-body Hamiltonian, and each fc-point will fol- 
low a circle in the Brilliouin zone. If the loop encircles 
the Dirac point, non-adiabatic charge pumping should 
take place, in which the wave function acquires a non- 
trivial AA phase. So the question is: can this be detected 
with a DC transport measurement? In order to answer 
this, we have formulated the Kubo formula for the DC 
response for systems in intense AC fields, which is ac- 
complished with the Floquet-matrix formalism by which 
we can solve the time-dependent dynamics of fc-points 
within a static approach. We shall show that a TKNN- 
like formula for the Hall conductivity is obtained, where 
the Berry curvature is now expressed in terms of Floquet 
states which depends on the AA phase. We apply this 
formula to a single Dirac band first, then a graphene. 
For graphene 0, 0], which is dominated by the chirality, 
we conclude that a photo- induced Hall current (despite 
the absence of uniform magnetic fields) should appear 
in graphene irradiated by circularly polarized light and 
attached to two electrodes, where the Hall current can 



exceed the longitudinal current in magnitude. 

Kubo formula in the presence of strong light fields - 
We first derive the Kubo formula for electric transport 
for systems in strong AC fields, where we concentrate 
on the one-body problem for simplicity. The AC electric 
field is introduced as a time-dependent gauge potential 
Aac(t), which satisfies A ac (t + T) = A ac (t) with T the 
periodicity (while the frequency is = 2ir/T), On top of 
this we introduce, as in the linear response, a weak gauge 
potential A(t) = Et that changes slowly to represent an 
infinitesimal DC electric field E, where we have set e = 
1, h= 1. Thus we have a time-dependent Hamiltonian, 

/dk 
^(k)h(k + A ac (t) + A(t))^(k), (1) 

where h(k) is the one-body Hamiltonian and \I> a 
state vector (which is multi-dimensional for multi- 
bands). To represent the states in AC fields, 
we can employ the Floquet operator (see [1, Q) 
H(fe,A ac (t),A(t)) = h(k + A ac {t) + A{t)) -idt, with 
which the time-dependent Schrodinger equation reads 
H{k,A ac (t),A(t))\^(k;t)) = 0. Since A(t) is infinites- 
imal and adiabatically changing, we take it as the adia- 
batic parameter, while the AC field is intense and rapidly 
oscillating. So, for each interval of time over which A 
may be considered to be constant, we can introduce the 
Floquet states (a time-analog of Bloch states) satisfy- 
ing the Floquet equation H(k, A ac (t), A)|$ Q (fe; A,t)) = 
e a (k;A)\& a (k; A,t)} with a periodicity \$ a (k;A,t + 
T)) = \$ a (k;A,t)) P3, where e a (k;A) is called the 
Floquet quasi-cnergy which is a sum of the dynamical 
phase and the AA phase (see eqn.JHJ) below), and a la- 
bels the eigenstate. The solution of the time-dependent 
Schrodinger equation for a fixed A can be expressed as 
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If we define an inner product 

averaged over a period by ((a\/3)} = ij, dt(a(t)\P(t)), 
the Floquet states form an orthonormal basis, i.e., 
(($ Q (fe; A)\$p(k; A))) = 5 a p. With these as a basis the 
solution to the time-dependent Schrodinger equation for 
the slow change in A is 

[*(fc;A(t),t)> = e -ifZdt'e a (k;A(t')) [|$ a (fc ; A(t), t))+ 
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up to first order in time derivatives, with the a-th Flo- 
quet state taken to be the initial state. One can readily 
derive this result with the two-time methodfllL [l2j. We 
can immediately notice that the geometrical phase ap- 
pears in a form A(t))|^ • ^ |$ Q (fc; A(i))» = 

£^ Q (<<Mfc;A(t))|M • f£|* (fc;A(f))»MA(t)) - 
e a (A(t))] _1 . Since the current operator is J = dh(k + 
A ac (i) + A)/9A, the above formula for the DC transport 
in an intense AC background field is rewritten as 



Cab(A ac ) 



dk v [fp(k) - f a (k)] 



^ (($ a (fc)ij b i^(fc)))((^(fc)ij a i$ a (fc))) 

e/3(fe) - e a (fc) + i?7 

where f a (k) is the non-equilibrium distribution (occu- 
pation fraction) of the a-th Floquet state, r\ a positive 
infinitesimal, and we have put the perturbation A = 
as in the linear-response theory. The essential difference 
from the conventional Kubo formula in the absence of 
AC fields is that the energy is replaced with the Floquet 
quasi-energy, and the inner product with a time averaged 
one. We note that similar expressions were obtained by 
Torres and Kunold in their study of microwave-assisted 
zero-resistance states |l3|. The Hall conductivity can be 



further simplified to a TKNN-like formula, 
( T xv (A^) = e 2 l £/ Q (fc)[V fc x A a (k)] z , (3) 

where Aa{k) = — i(($ Q (fe)| V^|$ a (fe)}}. We note that 
if we separate the Floquet index into a = (i,m) where 
i labels the original band and m is the Floquet index, 
then A. a {k) is independent of m. However, the occu- 
pation f a (k) depends on both indices. In equilibrium, 
/(i,m)(fe) = S m0 f F D(Ei{k)) holds where Ei(k) is the 
energy of iih state and /fd the Fermi-Dirac distribu- 
tion. In non-equilibrium, however, the distribution is 
non-universal, and depends on the detail of the system 
such as how the electrodes are attached, etc., so that 
a case-by-case study should be needed to determine the 
distribution and hence the DC transport of the system. 

Application to a Dirac band — Now we study the ef- 
fect of AC fields on the two-dimensional Dirac band. 
The Hamiltonian is H{t) = r z v[k x + Al c {t)}<j x + v{kv + 
Ay_ c (t)](jy, where t z = ±1 labels the chirality [K and K' 
points in graphene), v is the velocity (set to v = 1 here- 
after), and <Ji the Pauli matrices. The circularly polar- 
ized field is given as (A x c , A^.) = A(cosH,t,smflt) where 
A = F/fl with F being the field strength. Here we ne- 
glect the momentum of light since « « c, and consider 
direct transitions, which is a situation different from the 
Volkov solution fl4||. With Fourier transformed Floquet 
states |$(i)) = J2 m e~"™ 0t |u™} the Floquet equation be- 
comes 
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FIG. 1: (a) A trajectory of k + A ac (t) around a Dirac point 
in a circularly polarizated light field, (b) The Floquet quasi- 
energy (black curves) plotted against k x with k y — for 
F/Q = 0.2. The color coding represents the weight of the 
static (m = 0) component, (c) Density of states for vari- 
ous field strengths, (d) The photo-induced Berry curvature 
[V fc x Aa(fe)] for a = (l,m) for F/Q = 0.2. The series 
a — (2, m) have a similar behavior with an inverted sign. 
Frequency is Q/v = 1 throughout. 



where the Floquet 
i jf dtH(t)e^ m -^ nt has H r 
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agonal components, whereas the off-diagonal components 
depend on r z , i.e., H mm+1 = , H™" 1 ' 1 
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for t z = 1 and H mm+1 = (°%) , ii™"- 1 = (°- A ) for 
t z = — 1. This defines an eigenvalue problem for a block 
tri-diagonal matrix that can be solved numerically with 
a truncation at certain |m|. 

As an important effect of the AC field, gaps open at 
uj = integer xCl/2 in the quasi-energy band structure 
(Fig. 1(b)), reflecting a Dirac-band analog of the AC- 
Wannier-Stark ladder, as also seen in the density of states 

A(u) = -U dk E a Im^Sfe * Fig. [He). The gap at 
uj = ±0/2 is the largest, which is related to one-photon 
assisted transport discussed later. We note that this gap 
also opens in the linearly polarized case as studied in [la ] . 
More importantly, a new gap opens at k = 0, e = in 
the circularly polarized case. Indeed, one can show that 
the solution of the time-dependent Schrodinger equation 
at k = is 



|* a (fc = 0,t)) oc e" 
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where the Floquet state is labeled by a = (i,m) with 
2 = 1,2 representing the upper and lower branches of the 
Dirac band, m the Floquet index. The quasi-energy is 
e a = g i + mn with ^ = V ™ F +" , e 2 = 
The a 



■ mVL with i\ = v — 7," ~ r " , e 2 
- (1,-1), (2,0) bands are direct descendants 
of the original Dirac bands, and the dynamical gap 2k 
opening between them is n = 



2 
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The dynamical gap first grows quadratically with A, 
2k ~ 2A 2 /il, followed by an asymptote 2k ~ 2 A — O. 
An important property of the quasi-energy is that it is a 
sum of the dynamical phase and the AA phase, i.e., 

s a = (($ a \H(t)\$ a }) +1 % A /T, (6) 

where the AA phase is given by 

7 ^ A = T(($ a \id t \$ a )) = ±tt {[A(A/nf + I]- 1 ' 2 - l} ,(7) 

where ± refers to a — (1, — 1), (2, 0). In the adiabatic 
limit (f2 — > with a fixed A) it approaches to =F7r. We 
note that only the fc-points with |fe| < A acquires the AA 
phase, since otherwise the Dirac cone is not encircled. In 
the Berry curvature for the a = (2, 0) Floquet state (Fig. 
[ljd)) there is a conspicuous peak around k = 0, 

[Vfe x A a (k)] z ~ ±r z ^K{\k\ 2 + K 2 y 3 / 2 , (8) 

where ± corresponds to a = (1, to), (2, m). In this 
expression, two geometric quantities appear, where the 
Berry curvature comes from the perturbative treatment 
of the weak DC electric field whereas the AA phase 
emerges because of the time-periodic dynamics of k- 
points in intense AC fields. Due to the factor r z , the 
contribution from K and K' points in graphene will can- 
cel with each other if the distribution is identical between 
them. However, we shall see that, if we apply a static DC 



bias across the system, chirality (the valley symmetry) 
is degraded, which will lead to a non-trivial curvature, 
hence to a Hall current. 

Keldysh approach to photovoltaic transport in graphene 
— So we move on to a Keldysh Green's function analy- 
sis of transport properties in a graphenc irradiated by 
a circularly polarized light and attached to two elec- 
trodes. The system is described by an action, S — 

Jc ^ (^graphene ~\~ ^mix ^ ^electrodes) ; where /^graphene — 

Ylijtj c l(idt — tije %Aij ^)cj is the tight-binding model for 
graphenc with a hopping tij = —w for nearest neigh- 
bors, while £ mix = J2 kr V£ ix [{a r k Y c r + h.c] represents 
the coupling between the electrodes and graphenc, with 
the spin degrees of freedom ignored. The AC field is in- 
troduced by A%(t) = (r t - rj ) ■ A ac (t) with A ac {t) = 
(F/il) (cos fit, sin fit), where F = eaE is the normalized 
field strength (a: lattice const.). We assume that the 
electrodes are described by a fermion operator a r (r £ 
{L, R} labeling the left and right electrodes), for which a 
Fermi-Dirac distribution (a r V) = f r = [e^-M + l]- 1 
is assumed with electrode-dependent chemical potential 
fj, r . They are related to the DC bias V across the elec- 
trodes by /i£ = V/2, /i_R = —V/2. Imposing a periodic 
boundary condition in the direction (y) of the graphenc 
ribbon, the Keldysh Green's functions for each momen- 
tum k y become a matrix labeled by the site in the x- 
direction (i = 1, . . . , N) and by the Floquet index. The 
Green's functions satisfy 



Gk y M G k v H \ (iu) + nfl + iv)S m nSij - (H)™ n (k y ) 



L /2 -ir L (l-2/ L (w + mfl)) \ (iV R /2 -iT R (l - 2f R (u 

-iT L /2 J + ° iN ° mn I -iT R /2 



where G KRA are the Keldysh, retarded, and advanced 
Green's function, respectively, T r cx |V^ ix | 2 is the imagi- 
nary part of the self-energy due to the sample-electrode 
coupling. The effective Floquet Hamiltonian is defined 
by (H) mn (k y ) = i J T dte^ m -^ nt H(k y ]A ac (t)), where 
the indices i, j are suppressed. The current between 
sites i,j is determined from the lesser component G < by 



k-. m,n J V 



-i(m— n)Ot 



(9) 



k y mn 

Jij(t)(Gf) mn; ji(uj) + Jji(t)(G k )m»;tiH 



where Jij = 5H /5Aij is the current operator. In practice, 
we do not calculate the Keldysh component but use the 
Keldysh equation (c.f. [3) to relate G< with G A ,G R 
obtained by diagonalizing the Floquet Hamiltonian. 



In the obtained current distribution (Fig[2fa)) the po- 
larized light induces locally circulating currents in the 
absence of the bias V across the electrodes. There is no 
net current in the y-direction as it should. This current 
resembles the orbital magnetism which was predicted to 
arise when a perturbation induce a gap in the Dirac cone 
[I?], EH- Indeed, the circularly polarized light in the 
present case opens a gap around the band crossing (arrow 
in Fig[3]^b)), which has led to a similar DC effect. 

A striking finding here is that, when we switch on the 
bias voltage V, we have a photo-induced DC Hall cur- 
rent as well as a longitudinal (|| x) current (Fig. [21(b)). 
The Hall current is natrually inverted when the right cir- 
cularly polarized field A ac cx (cos fit, sin Qt) is changed 
into the left polarization cx (cos fit, — sin fit), or the 
bias voltage is inverted. The / — ^-characteristics is 
shown in Fig. [3] (c) for J x and J y , the averaged cur- 
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FIG. 2: DC current distribution in an armchair graphene 
ribbon attached to electrodes subject to a circularly polarized 
light for a finite AC field F = 0.025u), Q = 0.3u; with no DC 
bias (a) and with a finite bias V = 0.005w (b). (c) / — 
V characteristics of the longitudinal J x (black) and the Hall 
current J y (blue) for various values of F. (d) DC conductance 
G — J/V plotted against field strength F for a fixed bias V = 
0.005w. System size in the x-direction N = 34 throughout. 

rent in x and y directions, respectively. The photo- 
induced net Hall current grows linearly with the bias V, 
but saturates and then decrease when V becomes large. 
Figured) depicts the dependence of the conductance, 
Gxx = Jx/V, G xy = Jy/V on the intensity F of the 
circularly polarized light for a fixed bias V . The conduc- 
tance, unlike the conductivity discussed in the first part 
of the paper, is less universal and depends on the contact 
r r (here we set Tl = Tr = Q.2w) etc, but we expect that 
the effect will be qualitatively robust. As we increase F, 
the longitudinal G xx first decreases, and increase again 
(F > 0.03w). The decrease can be explained by the gap 
opening at the Dirac points, while the increase is due to 



photo-assisted transport. We note that similar features 
have been experimentally observed in microwave irradi- 
ated carbon nanotubes The Hall conductance G xy , 
on the other hand, initially grows quadratically with F 
and then increases linearly, which is a dependence simi- 
lar to the gap k in the Dirac cone (eqn.©), which may 
indicate J y ~ kV . We note that a similar expression was 
obtained in the case where the chirality is broken in a 
static manner (20j. 

To summarize, we have found that a combined effect 
of an intense AC field of a circularly polarized light and 
a (weak) DC bias can produce a photo-voltaic dc Hall 
current in graphene, despite the absence of a uniform 
magnetic field. The typical intensity of laser conceived 
here, F ~ 0.001k;, corresponds to E ~ 10 7 V/m for pho- 




FIG. 3: Floquet quasi-energy diagram for an armchair [2lT] 
ribbon in an AC field for F = 0.005w(a), 0.02w(b) with fre- 
quency Q, = 0.3™. The color coding represents the static 
weight |(u° K)| 2 as in Fig. 



ton energy fl ~ leV, w = 2.7eV, a = 2.6A, which should 
be within the experimental feasibility. Inclusion of dissi- 
pation etc will be an interesting future problem. 
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